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THREE-DIMENSIONAL  ANALYSIS  OF  ACOUSTIC  SCATTTERING 
FROM  A  COATED  CYLINDRICAL  SHELL 


INTRODUCTION 

In  this  study,  a  theoretical  model  is  developed  for  evaluating  the  farfield  pressure  scattered 
from  an  infinitely  long,  coated  cylindrical  shell  by  a  plane  acoustic  wave  incident  at  an  angle 
with  the  normal  to  the  cylindrical  shell  axis  (a  three-dimensional  problem).  The  model  is  a  two- 
layer-structure  comprising  an  outer  layer  of  microvoided  elastomer  (coating)  that  is  perfectly 
bonded  to  a  cylindrical  shell.  The  coating  consists  of  an  acoustically  soft  material  that  is 
designed  for  reducing  flexural  wave  noise. 


BACKGROUND 

The  acoustic  reflection  from  a  solid  cylinder  in  water  is  described  by  Morse1  and  Bowen 
et  al.  The  scattering  of  sound  by  hollow  elastic  cylinders  in  water  has  been  investigated  by 
Gaunaurd3  and  Flax  and  Neubauer.4  In  these  cases,  the  incident  sound  wave  was  taken  as  an 
infinite  plane  wave  normally  incident  to  the  cylindrical  axis  (a  two-dimensional  problem).  This 
wave  excites  the  radial  and  circumferential  modes.  However,  a  plane  wave  incident  at  an  angle 
with  the  normal  to  the  cylindrical  axis  excites  modes  of  vibration  of  the  cylindrical  shell  having 
radial,  circumferential,  and  axial  dependence  (a  three-dimensional  problem). 

SPECIFIC  OBJECTIVE 

This  report  describes  the  three-dimensional  analysis  for  the  two-layer  coated  cylindrical 
shell.  An  earlier  two-dimensional  problem5  has  been  extended  to  the  three-dimensional  case  to 
investigate  the  scattering  of  the  plane  acoustic  wave  by  the  coated  cylindrical  shell.  Ko  and 
Sandman6  presented  a  portion  of  the  two-dimensional  results.  The  formulation  of  the  problem  is 
based  on  the  theory  of  elasticity,  acoustic  wave  equations,  and  pertinent  boundary  conditions. 
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The  major  accomplishment  of  this  research  is  the  numerical  calculation  of  directivity  patterns  for 
the  coated  cylindrical  shell  using  the  three-dimensional  analysis. 
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THEORETICAL  ANALYSIS  OF  ELASTIC  WAVES 
IN  CYLINDRICAL  COORDINATES 


The  geometry  of  the  theoretical  model  in  this  study  is  depicted  in  figure  1.  The  outer 
surface  of  the  composite  structure,  which  consists  of  the  coating  and  the  cylindrical  shell,  is 
contact  with  water.  The  core  (cavity)  of  the  structure  contains  air. 


Figure  1.  Geometry  of  Theoretical  Model 

It  is  essential  to  study  waves  propagating  in  elastic  media.  In  this  section,  equations 
necessary  for  the  formulation  of  the  present  problem  are  developed.  The  vector  differential 
equation  that  governs  the  small  elastic  motion  in  the  elastic  medium  is  written  as7,8 


o 
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(1) 


/iV2K  +  (l  +  /i)V(V.«)  =  p§, 

01 

where  X  and  p  are  the  Lame  constants,  u  is  the  displacement  vector,  p  is  the  material  density,  V 

2 

is  the  gradient  operator,  t  is  the  time,  and  V  is  the  Laplacian  operator.  The  solutions  of 
equation  (1)  are  a  combination  of  a  vector  potential  ^and  a  scalar  potential  function  <f>,  so  that 


«  =  V^  +  (Vxf), 


(2) 


which  is  true  provided  that  <f)  and  ^are  the  solutions  of  the  elastic  wave  equations 


’•'■if' 


(3) 


Y72  1  5V 


(4) 


and 


V-y/  =  F, 


(5) 


where  F  is  an  arbitrary  function  of  the  spatial  coordinate  and  time.9’10  The  complex  dilatational 
(compressional)  and  shear  (transverse)  wave  speeds  are  given  by 
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X  + 


-|l/2 
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and 
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If  equation  (3)  is  written  in  cylindrical  coordinates,  then 


d2$  \  d<j>  1  d2<f>  d2(f>  _  1  d2<j> 
dr 2  r  dr  r2  d62  dz2  c2d  dt 2  ’ 


where  r,  0,  and  z  are  the  spatial  coordinates  in  the  radial,  circumferential,  and  axial  directions, 
respectively.  If  equation  (2)  is  used,  then  the  radial  displacement  is  written  as 


»,  =  (V*  +  Vx*0,=|£+ 

dr 


1  dyrt  dy^g  N 
r  d0  dz  j 


(9) 


where  the  radial  displacement  ur  is  shown  in  figure  2.  Note  that  other  components  of  the 
displacement  vector  and  related  stress  tensors  are  also  shown. 

Similarly,  the  circumferential  and  axial  displacements  are  written  as 


r  od  \  dz  or 


(10) 


Figure  2.  Displacement  Vector  Components  and  Stress 
Tensors  in  Cylindrical  Coordinates 
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and 


u2  =  (V*  +  Vx*r)f 

(11) 

=  M  +  fe  +  ^£  1 

3z  r  dr  r  36  / 


Equation  (4)  can  be  decomposed  in  the  radial,  circumferential,  and  axial  directions.  In  the  radial 
direction,  one  obtains 


2  3^  _  i  ay, 

r2  30  c2  3/2 


r&¥r  |  1  !  1  9Vr  ,  ^Vr  ' 

k  3r2  r  a-  r2  302  3z2  , 


2  dyg  _  1  3>r 
r2  30  c2s  dt 2 


(12) 


In  the  circumferential  direction,  one  obtains 


2  dy/r  _  1  3>e 
r 2  30  c2  dt2 


ff gVg  ,1^8,  1  g¥g  gW|  jgg.1  2  3^r  _  1  3>g  (13) 

3r2  r  dr  r2  302  3z2  J  r2  J  r2  30  dt 2 

In  the  axial  direction,  one  obtains 


vy2  = 


_  i  aV. 


3/ 2 


3V,  1  3^.  1  3V,  3V,  1  32^ 

"a^rlT+r5’  302  +  ~dzr~~fs~dF 


(14) 
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The  solution  of  equation  (8)  is  written  as 


where  k,  =  k0  sin^?  is  the  wave  number  in  the  z-direction,  k0  =  co  lc0  is  the  acoustic  wave 
number,  co  is  the  frequency  in  radians/s,  and  c0  is  the  sound  speed  in  water.  Substituting 
equation  (15)  into  equation  (8)  gives 


^d2<&  loT<I>  1  d20^ 

0— t-  +  — — 0  +— — 


Q)‘ 


dr  r  dr 


r  ddz 


■k2<$>Q)  =  -—7-00. 


Dividing  equation  (16)  by  0,  one  obtains 


d  $>  1  d<t>  2 . 

— r  + - +  a  <■ I>  = 

dr  r  dr 


rd2®^ 

d92 


0 


<6, 


where  a  ~(k2  -  k2)'2,  and  kd 


co  I cd  .  Dividing  equation  (17)  by  O/r2  gives 


d2< D 
dr 


■  +  r 


d® 

dr 

<D 


2  2 

+  a  r 


d 20 

del 

0 


Let 


d20 

*E=k2e 

©  * 


Then,  multiplying  equation  (IS)  by  <t>/r2y  one  obtains 


0  =  0. 


d2<$>  1  d<&  (  2 

i/r  r  dr  ^  r  J 

(20) 

Equation  (19)  can  be  rewritten  as 


d2<d 

~d¥ 


+  k2g®  =  0. 


(21) 


The  solution  of  equation  (21)  is  written  as 
0(0)  =  A  cos (ke0)  +  B sin(ke0) . 
Single-valueness  requirements  on  0(0  give 
ke=n  (an  integer). 


Then,  equation  (20)  becomes 


d 2<D  \d$  (  2 

— r+ - +  a 

dr  r  dr  y 


0  =  0. 


(22) 


(23) 


(24) 


Note  that  equation  (24)  is  the  Bessel  equation  of  order  n. 

Equation  (15)  is  now  written  as 

<}>{r,0,z,(p,f)  =  <t>(r)cos,n@  (25) 

Then,  the  solution  of  equation  (8)  is  written  as 
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Mr,f),z,(p,t)  =  {A^Jarj  +  B/Jar^cosnO  e 


(26) 


-i(k.z-(ot) 


and  J nicer)  and  Y„(ar)  are  the  first  and  second  kinds  of  the  Bessel  function  of  order  n,  with  the 
argument  ar,  respectively. 

Similarly,  the  solutions  of  equations  (12)  and  (13)  are  assumed: 


yfr(r,0,z,<p,t)  =  '¥r(r)smnOe-i(k*z-0,t) 


(27) 


and 


V6ir,e,z,<p,t )  =  %(r)cos  nOe**-* 


(28) 


Substituting  equations  (27)  and  (28)  into  equation  (12)  gives 


d2Vr  1  d¥r 

- —  H - - 

dr 2  r  dr 


+  ^(- n2%  -Vr  +  2n%)+02% 


=  0, 


(29) 


where  /?  =  (k2  -  k2  Jn  and  ks=co/cs. 


Similarly,  substituting  equations  (27)  and  (28)  into  equation  (13)  gives 


d2Vff  1  d¥a  1 


+  +2nVr)+p2%=0. 


dr  r  dr  r 


(30) 


Subtracting  equation  (30)  from  equation  (29)  gives 


'  d2  1  d 
<{  - ~  + - + 

dr  r  dr 

L  >*  JJ 

•('P,-%)=0. 


(31) 
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Adding  equation  (30)  to  equation  (29)  gives 


\  d2  1  d 
[dr  r  dr 


(”-l)2 

2 


(X  +  X)=0. 


The  solution  of  equation  (31)  is  written  as 


X  -  X  =  2 (fir) + 2B,r„,  (fir) 


Similarly,  the  solution  of  equation  (32)  is  written  as 


(32) 


(33) 


V„+%  =24j,.,(fir)  +  2B'2r,-,m.  (34) 

The  property  of  the  gauge  invariance9, 10,11  can  now  be  utilized  to  eliminate  two  of  the  integration 
constants.  It  may  be  shown  that  any  one  of  three  potentials,  <£>,  (Yr  -  Ye),  and 
(Yr  +  Y*),  can  be  set  equal  to  zero,  without  loss  of  the  generality  of  solution.  If  (Yr  +  Ye)  =  0, 
one  obtains 


x = -x , 


(35) 


which  also  means  that  A2  =  B2  =  0 . 

Substituting  equation  (35)  into  equation  (33)  gives 

X(.r)=^JM(fir)^B,Y„t(fir).  (36) 

Note  that  Ye  (r)  =  -Yr  (r) ;  i. e., 

X(r)  =  -AJ„,m-BX.„(firy  (37) 
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If  equation  (36)  is  substituted  into  equation  (27),  then  the  radial  component  of  the  vector 
potential  is  written  as 


V.  (T,  e,  Z,  <P,  t )  =  fa  J.tl  (fir)  +  BJ„,  (fir)}  sin 
Similarly,  the  circumferential  component  of  the  vector  potential  is  written  as 
MrAz,9A  =  {-[AJ,.,(fir)+BJ,tfifir)}}casn8e-‘^--'\ 

Finally,  the  axial  component  of  the  vector  potential  is  assumed  to  be 
y/2  (r,  9, z,  <p,  t)  =  %  0 r ) sin  n6  e-i(k^wt) . 

Substituting  equation  (40)  into  equation  (14)  gives 

X(r)  =  A3Jnm+BiYn(M- 

Combining  equation  (41)  with  equation  (40)  gives 
^  (r,0,  z,^,0  =  {V„  G^)  + 

If  equations  (26),  (39),  and  (42)  are  substituted  into  equation  (9),  the  radial  displacement  is 


written  as 


dr  r  ee  dz 


-4 


a-J,(ar) 


+  B1 


ccYn(ar ) 


+ 


-ik-Jn+i(Pr) 


“h  Bn 


+  A 

+  B 3 

~Yn(Pr)]> 

__r 

J  J 

(43) 


xcos  nOe-^™0. 


If  equations  (26),  (38),  and  (42)  are  substituted  into  equation  (10),  the  circumferential 
displacement  is  written  as 


eK  ’  ’  J  r  dG  dz  dr 


A 


n 


•/„(«') 


+  Bl 


n 


~~Yn(ar) 


+A 


ik:Jn+l(Pr) 


+  B~ 


-ikJn^Pr) 


r 

> 

+  A 

-pj'ntfr) 

+  B3 

-PY'SPv)  > 

j 

(44) 


xsinw^e-'^0. 

If  equations  (26),  (38),  and  (39)  are  substituted  into  equation  (1 1),  the  axial  displacement  is 
written  as 
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z  5z  r  dr  r  d6 


A 


ikzJ „{ar) 


+  A 


-ik  ^(ar) 


4*  A 


L  \  r  J 


+  Bn 


\  r  J 


xcos  n$e<k*-*'\ 


(45) 


The  radial,  circumferential,  and  axial  stress  components  are  written  as 


r  =  A 


5m_  1  du,  ur  du , 

— c  + - S.  +  -L  +  .  * 


dr  r  dO  r  dz 


^  o  dii  r 

+  2/i— L 


Sr 


=  (A  +  2/of“l  +  [(A  +  2//)-2//]  fl^L  +  ^  +  f^O 

\3rJ  r  az  J 


(46) 


f  1  a«r  due  Un  ^ 

r,  =  u - r-  +  — 6-  -  -J- 

\r  d9  dr  r  j 


(47) 


and 


*rz=B 


du.  dur  N 
dr  dz  j 


(48) 
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Substituting  equations  (43),  (44),  and  (45)  into  equation  (46),  one  obtains  the  radial  stress 


rjr,  9,z,<p,t) 


Pc] 


(duA 

v  dr  ) 


\r  o9  r  dz ) 


AApcd-a  J „  (ar) 


p(c2d-lc]) 


'i,V 

<r) 


+  k]  \Jn(ar)  +  jjl(ar) 


hi 


+  Bx\pc]-ct2Y"n(ar) 


2  cfl 

—  < 

+  k]\Y„(ar)  +  -rl(.ar)  ■ 

y) 

J  r  J 

+  A2\pc][-ikJjUf3r)) 

+  pfe-2  c*)[ik,pfUpr)^ 
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+  B. 


+ 


p{cl-2c])\ikJY'M(Jr)\ 


+  A,  [pc]  {-  4  J,  (fir)  +  IfiJl  (fir)' 
{  v  r  r  J 


+  P 


te-2 *i) 


-^-f.(fir)  +  ^J„(fir) 
r  r 


+B,\p4~Y.  (fir)  +  -fit,  (fir) 

l  V  r  r  J 


p(c]  -2c;) 


n±Y',(fir)^Y,(fir) 
r  r 


:cos  nde-Kk‘*-at\ 


Substituting  equation  (43)  and  (44)  into  equation  (47),  one  obtains  the  circumferential  stress 


i:re(r>0>z><P’t)  =  Pc. 


1  dur  dUg  Ug 

r  d9  dr  r 


A\pc. 


2 n  2 na  < 

—JAar) - Jn(ar) 


+  BA  pc] 


2 n  2na  •  ' 

~^Y„(ar) — —  Yn(ar) 


+  M  P°l 


i(n  + 1  )k2 


(49) 
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+ BA  p&s 


i(n  + 1) 


kzYn+l{Pr)-ikJYn^fir) 


+  A.  i 


pc ; 


f rp2 

UJ 


r  J 


, 

f  wV 

< - - 

*<■> 

r, (fir) -  fi%(fir)  S-Ufir) 

r  J 


xsin  ndeKk^at) . 

Substituting  equations  (43)  and  (45)  into  equation  (48),  one  obtains  the  axial  stress 

zrz(r,e,z,<p,t)  =  pcs\ -zr-  +  —L- 

lz  J 


pc2 

Sr 

r 

T 

4 

k 

L 

-lilc.aJ  n(ar) 


+  B, 


PCs 


-  2ikzaYn  ( ar ) 


+  4i  pc] 


+  B2  \pC] 


+  A{  Pc: 


+  BA  pel 


'(g+j)  tk 

2  Kz 

\  r  j 


'(«  +  !) 
^  r2 


-kt 


(fir)  ~  fiJi ,  (AO  -  A’C.  (AO 

r 

n„  (^)  -  AC  (AO  -  A’ Cl  (AO 


ink.  T  /  n  . 

- ^„(/fr) 

r 


ink. 


y,(fir) 


(50) 


(51) 


xcosnfe-''1''-*"’. 


16 


The  divergence  of  y/is  an  arbitrary  function:9,10 


V  \j/  = 


dr 


Vr  !  1  dVe 
r  r  36 


+ 


Q¥= 

dz 


(52) 


In  the  outside  fluid  medium  (water),  the  acoustic  pressure  is  given  by  the  sum  of  the 
incident  wave  and  the  scattered  wave  pressures  as  follows: 


Po  (r,  6,z,<p,t)  =  pt  (r,  6,z,(p,t)+  ps  (r,  6,  z,  <p,t),  (53) 

where  p0(r,6,z,(p,t)  is  the  total  wave  pressure,  Pi(r,6,z,(p,t )  is  the  incident  wave  pressure,  and 
ps  (r,  6 ,  z,<p,  t )  is  the  scattered  wave  pressure.  When  a  plane  wave  is  incident  at  an  angle  (p  with 
the  normal  to  the  cylindrical  axis,  the  total  pressure  is  written  as 


=  Pi e',(k^at) £ (-/)X |  J„  (Vcosp)  +  Aln)H™{k0r cos (p^cosnO,  (54) 

where  P,  is  the  amplitude  of  the  incident  wave,  A^n)  is  the  unknown  coefficient  to  be  determined 

for  order  n,  and  sn  is  the  Neumann  constant  (sn  =  1  for  n  =  0  andf„  =  2  for  n>  1).  Note  that  the 

scattered  outgoing  wave  that  satisfies  the  radiation  condition  is  given  by  the  Hankel  function  of 
the  second  kind: 

H^\k0rcos(p)  =  Jn(kQrcos(p)  -  iY„  (k0rcos<p).  (55) 

Inside  the  coated  cylindrical  shell  (air),  the  pressure  field  is  written  as 

p,(r,e,z,<p,t)  =  P,  cos nff ,  (56) 

n-  0 
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where  p(j,  9,z,  <p,f)  is  the  pressure  field  in  the  core  of  the  cylindrical  structure,  B(0n)  is  the 
unknown  coefficient  to  be  determined  for  order  n,  ki=a>lci,ci  is  the  sound  speed  in  the  core, 
and  cos#>3  =  [1  -  (k0  sin q>/k3)2 ]1/2 . 

To  ensure  consistency  with  equations  (54)  and  (56),  the  scalar  potential  <j)(r,  9,z,  (p,t)  and 
the  three  components  of  the  vector  potential  y/(r,  9,  z,<p,t)  may  be  rewritten  as 

(_ox  (ar)+B<’>y„  ( ar )}  cosnS,  (57) 

n= 0 

v,  0 - ,  o,z,<p,t)=pi  (-0"  £,  {4 *v„,  m + wo}  sin  „e,  (58) 

n=0 

%  (r,e,z,vA  =  Ple^-“'%-lye,t-4’>J,,l(J)r)-B™r„,(f)r)}  cosnff,  (59) 

n= 0 

and 

V,  (.r,e,z,9,t)  =  7>e-,<t-'-'>i;(-/)V.{4”>J„WO+Ss<">i;  (fir)}  *■»«  (60) 

*= 0 

If  equations  (57)  through  (60)  are  used,  the  radial  displacement  ur,  the  circumferential 
displacement  ue,  the  axial  displacement  uz,  the  radial  stress  xn,  the  circumferential  stress  rr6, 

oo 

and  the  axial  stress  zn,  are  preceded  by  the  expression,  (-/)”£■„ .  These  six  quantities 

n~  0 

include  the  six  unknown  coefficients:  A™ ,B[n) ,A^n) ,B^ ,A^n) ,  and  B^n) .  These  coefficients 
should  be  determined  for  each  order  n  by  using  the  pertinent  boundary  conditions. 
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FORMULATION  OF  THE  PROBLEM 


In  this  section,  the  unknown  coefficients  are  determined  by  using  the  pertinent  boundary 
conditions.  The  boundary  conditions  to  be  satisfied  at  the  interface  between  the  outer  fluid  and 
the  outer  surface  of  the  coating  (denoted  by  layer  1)  of  the  cylindrical  shell  structure;  i.e.,  r  =  R0, 
are  written  as 


]r=R0  ’ 


(61) 


f0V?1 

1 

[>ri 

1 - 

9> 

to 

L _ 

r=R„ 

dr 

(62) 


(63) 


and 


[rrej]r=s0- 0,  (64) 

where  [r^]r=i?o  is  the  normal  stress  of  order  n  at  the  surface  r  =  R0 ,  and  the  subscript  1  refers  to 

the  layer  1  (coating).  Similar  notations  are  used  for  the  total  pressure,  the  radial  displacement, 
and  the  shear  stresses.  Note  that  p0  is  the  density  of  the  outer  fluid.  The  boundary  conditions  to 

be  satisfied  at  the  interface  between  the  cylindrical  shell  and  the  coating,  i.e.,  r  =  Rv  are  written 
as 


k”i]r -arteL*. 


(65) 
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k’,L,= 

teL,, 

(66) 

kiU= 1 

;  rz^l\  r=i?j  > 

(67) 

k?L,=  l 

W(W)] 

.  r,2  Jr^  ? 

(68) 

["»”’]»*  =  1 

>)] 

“0,2  J  r=R1  ? 

(69) 

and 

k,’U=k?U-  po) 

The  boundary  conditions  to  be  satisfied  at  the  interface  between  the  cylindrical  shell  and  the  core 
fluid,  i.e.,  r=R.2,  are  written  as 


kZLR=  0,  (71) 

te]r^=0,  (72) 


'd2u(rf 

1 

r*f”i 

cS 

_ 1 

r=J?2 

.  &  . 

(73) 


and 


]  r=J?2  —  [  Pi  }]r=R2  >  (74) 

where  pin)  is  the  pressure  field  in  the  core,  and  p%  is  the  density.  Substituting  equations  (49) 
and  (54)  into  equation  (61)  gives  for  order  n 
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(75) 


a,(r)4o")  +«i(2^u  +4”Xi)+fli(4)4”) +«,W+«£)4g)+flft)^;) = b\n\ 


where  the  superscript  (n)  denotes  the  order  n,  and  the  coefficients  ... ,  and  £,(n)  are 

given  in  appendix  A.  Note  that  A[f  is  the  unknown  coefficient  A  i  of  order  n  for  the  layer  1 

(coating).  Similar  notations  are  used  for  other  unknown  coefficients.  Substituting  equations  (43) 
and  (54)  into  equation  (62)  gives  for  order  n 


+  a™ 4$  +  a^B[f  +  +  agBjg  +  agA™  +  a%>B$  =  i<"> . 


(76) 


Substituting  equation  (50)  into  equation  (63)  gives  for  order  n 


ag’4?  +<5,?  +aj?4J  +ag>4j>  =  0, 


(77) 


Substituting  equation  (51)  into  equation  (64)  gives  for  order  n 

a$>4?+a <;>£!;>  +a£>4")  +a«4’>  +a£>jj«  =  0. 


(78) 


Substituting  equation  (49)  into  equation  (65)  gives  for  order  n 

<4?  +«£’*&’  +<4?  +<4?  +<4?  +4"’4?  +<4? 
+<4?  +4S45  +4S«g +4S4?  +454?  =o. 


(79) 


Substituting  equation  (50)  into  equation  (66)  gives  for  order  n 

<4?  +  454?  +*£’4?  +454?  +454?  +o«B<->  +a«4? 
+454?  +a«45  +454?  +a«4”>  +a«B«  =  0. 


Substituting  equation  (51)  into  equation  (67)  gives  for  order  n 


(80) 
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4545  +  4545  +  454? + <45  +  454? + <45  +  <4? 
+454?  +45,45  +45,45  +45  4?  +4?, 4?  =  0. 

Substituting  equation  (43)  into  equation  (68)  gives  for  order  n 


4545 +<45  +<45+4545  +<45  +4545+454? 
+<4? +<.4?  +4545  +4545  +4345  =  o. 


Substituting  equation  (44)  into  equation  (69)  gives  for  order  n 


+4545  +<4?  +<45  +<4?  +434?  +<4,5 
+<45  +45.45 + 45,45  +4545 +4545 = 


Substituting  equation  (45)  into  equation  (70)  gives  for  order  n 


„0)  a(»)  .  „(«)  dO 


(«)  I  nw 
'  ^10,4 


*10,5^2,1 


^w10, 8^1,2  ^  “10,9^1,2  ^^10,10^2,2  ^w10,11^2,2  ” 


Substituting  equation  (50)  into  equation  (71)  gives  for  order  n 


✓,00  /f(»)  ./,(«) 
^11,8^1,2  '  ^11,9-°!,: 


(«)  l/y(»)  jO) 
2  tw1U0^2,2 


+45,45  +45.45  +45,45  =  o. 


Substituting  equation  (51)  into  equation  (72)  gives  for  order/? 


+ 42,45 + 42,45 + 42,45  =  ». 


Substituting  equations  (43)  and  (56)  into  equation  (73)  gives  for  order  n 


Substituting  equations  (49)  and  (56)  into  equation  (74)  gives  for  order  n 


4?+43-bS’+<!.4? 


+<OS  +<C4? +4S.4? +<1.4”  =  o. 


(88) 


Using  equations  (75)  through  (88),  one  may  obtain  a  system  of  linear  algebraic  equations  to  be 
solved  for  the  unknown  coefficients  for  order  n:  A(0n),  Af”\  B[f,  B^,  A^,  B$,  A$ , 

B$,  A^,  B^,  A^,  B^,  and  B[n)  (shown  in  appendix  B).  After  the  coefficient  A(0n)  has  been 
obtained  by  solving  the  algebraic  equations  shown  in  appendix  B,  the  pressure  field  scattered 
from  the  coated  cylindrical  shell  at  distance  r  is  written  as 


Ps  (r,  0,  z, <p, t)  =  (-/)"  (k0r  cosp)  cos n9 . 


(89) 


n= 0 


Further,  the  farfield  pressure  can  be  obtained  by  using  the  asymptotic  form  of  the  Hankel 
function.  Its  asymptotic  form  (r-»oo)  is  given  by 


,1/2 


H(n2\k0rcos<p) 


xk^rcoscp  J 


exp 


f ,  nn  7r^ 

k0rcos<p~— —  — 

V  2  4yj 


(90) 


If  equation  (90)  is  combined  with  equation  (89),  the  farfield  pressure  is  written  as 


PM<P)  =  PiQx P 


-/ 


7t 


k0rcos<p-k  z - 

4y 


,1/2 


7rk0rcos(p  J 


2' s^cosnO . 


n-0 


(91) 


23/(24  blank) 


NUMERICAL  CALCULATIONS  AND  DISCUSSION 


The  numerical  calculations  for  the  scattered  field  from  the  coated  cylindrical  shell  have 
been  made  using  equation  (91).  The  normalized  farfield  pressure  level  as  a  function  of  9  for  a 
given  (p  is  expressed  as 


x¥(9,<p)  =  20  log. 


0(0, P) 

_®(0,P)  max 


dB, 


(92) 


where 


f  2  'l 

1/2 

+oo 

®(0,p)=/i 

^snAgn)  cos»<? 

yxk0r  cos  <p  y 

«= 0 

(93) 


^{9,<p)  is  the  directivity  pattern  for  the  coated  cylindrical  shell,  A(0n)  is  obtained  from  the 
solution  of  the  equation  shown  in  appendix  B,  and  (i>(9,^>)max  is  the  maximum  value  of  <t>(9,<p) 
calculated  for  9=0-  360° using  equation  (93).  The  directivity  pattern  for  a  rigid  cylinder  can  be 
obtained  from  equation  (92)  provided  that  the  coefficient  A[n)  is  given  by 


/(M  _  UKR) 

^0  _  I  5 

(M) 


(94) 


where  R  is  the  radius  of  a  rigid  cylinder.3’4 

In  calculating  the  directivity  patterns  using  equation  (92),  it  is  necessary  to  use  the 
properties  of  the  coating  and  cylindrical  shell,  as  well  as  their  dimensions.  The  major  parameters 
required  in  the  present  study  are  the  material  density,  the  dilatational  (compressional)  and  shear 
(transverse)  wave  speeds  in  the  material,  and  the  densities  and  acoustic  wave  speeds  of  the  fluid 
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media  in  contact  with  the  outer  and  inner  surfaces  of  the  cylindrical  shell  structure.  Note  that 
dilatational  and  shear  wave  speeds  are  given  by  equations  (6)  and  (7).  The  Lame  constants  X  and 
/u  shown  in  both  equations  can  be  expressed  in  terms  of  elastic  constants  as  follows: 


A=£<r/[(l  +  o-)(l-2a)]  and  /*=£/[2(l  +  cr)],  (95) 

where  E  is  Young’s  modulus  and  cris  Poisson’s  ratio.  The  elastomer  coating  (rubber-like 
material)  has  the  Young’s  modulus  ranging  on  the  order  of  10  to  10  dynes/cm  as  a  function  of 
frequency.  The  Poisson’s  ratio  crfor  the  rubber  material  normally  approaches  0.5.  Rubber  is  a 
lossy  material  that  has  an  imaginary  part  of  the  Young’s  modulus.  The  complex  (dynamic) 
Young’s  modulus  is  given  by  E  =  Er+iEi,  where  Er  and  Ei  are  the  real  and  imaginary  parts  of 

the  Young’s  modulus,  respectively.  Actual  values  of  Er  and  Ei  for  a  given  material  can  be 
obtained  from  measurements.  The  ratio  of  Et  to  Er  for  rubber  ranges  from  approximately  0. 1  to 
1 .0,  depending  on  the  frequency.  As  can  be  seen  in  equation  (95),  X  cannot  be  defined  for 
rubber-like  materials  because  X  is  infinite  when  cr=  0.5.  Therefore,  the  dilatational  wave  speed 
(cj)  cannot  be  obtained  from  equation  (6),  but  the  shear  wave  speed  (cs)  can  be  obtained  from 
equation  (7).  An  alternative  expression  for  the  dilatational  wave  speed  in  rubber  is  written  as 

cd=[(B  +  4ju/3)/p]m,  (96) 

where  B  is  the  complex  bulk  modulus  of  material.  Note  that  |f?| » \ju\.  The  loss  factor 

associated  with  the  dilatational  wave  is  normally  very  small  and  is  weakly  dependent  on 
frequency.  The  complex  dilatational  and  shear  wave  speeds  for  the  layer  1  (coating)  are  now 
written  as  follows: 

cdX=cdm{\+Kj'2  (97) 


and 
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c,i  -^oiO+'Ci)  > 


(98) 


where  cd0l  and  ci01  are  the  dilatational  and  shear  wave  speeds  in  the  coating,  and  and  C,s\  are 

the  loss  factors  associated  with  dilatational  and  shear  waves,  respectively.  The  cylindrical  shell 
used  in  the  present  study  is  steel.  The  Lame  constants  X  and  p  for  the  cylindrical  shell  can  be 
expressed  in  terms  ofE  and  eras  shown  in  equation  (95).  Normally,  the  Young’s  modulus  of 
steel  is  real  because  it  is  not  considered  to  be  lossy  material.  However,  all  real  structures  possess 
some  inherent  structural  damping.  Thus,  a  frequency-independent  small  loss  factor  can  be 
assigned  in  the  numerical  calculation.  Then,  the  complex  dilatational  and  shear  wave  speeds  in 
the  cylindrical  shell  can  be  directly  obtained  by  using  equations  (6)  and  (7)  in  terms  of  the 
complex  Young’s  modulus  and  Poisson’s  ratio  and  are  written  as  follows: 

Cd2  ~  Cd02  0  +  d2 )  >  (99) 


and 


Cs2  ~  CsOl0-  +/£j2)  >  (190) 

where  cd02  and  cs02  are  the  dilatational  and  shear  wave  speeds  in  the  layer  2  (shell),  and  £d2  and 
CS2  are  the  loss  factors  associated  with  dilatational  and  shear  waves,  respectively.  The  baseline 
data  used  in  the  calculation  of  the  directivity  patterns  are  as  follows: 


R2  (inner  radius  of  cylindrical  shell) 

25.4  cm 

h2  (shell  thickness) 

5.08  cm 

\  (coating  thickness) 

5.08  cm 

p0  (water  density) 

1.0  g/cm3 

c0  (sound  speed  in  water) 

150  000  cm/s 

p3  (air  density) 

0.00121  g/cm' 

c$  (sound  speed  in  air) 

34  000  cm/s 
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p,  (coating  density) 

0.6  g/cm 3 

csQ]  (shear  wave  speed  in  the  coating) 

5000  cm/s 

cd01  (dilatational  wave  speed  in  the  coating) 

20000  cm/s 

£sl  (shear  loss  factor  of  the  coating) 

0.3 

£d]  (dilatational  loss  factor  of  the  coating) 

0.03 

p2  (shell  density) 

7.8  g/cm3 

E2  (Young’s  modulus  of  shell) 

19.5  x  10n  dyn/cm2 

a 2  (Poisson’s  ratio  of  shell) 

0.3 

£s. 2  (shear  loss  factor  of  shell) 

0.01 

C,d2  (dilatational  loss  factor  of  shell) 

0.001 

Figure  3  shows  a  comparison  between  the  directivity  patterns  calculated  at /=  5000  Hz 
using  the  two-  and  three-dimensional  analyses.  The  solid  line  shows  the  directivity  pattern 
calculated  using  the  two-dimensional  model  and  that  calculated  using  the  three-dimensional 
model  with  <p=  0.  As  shown  in  this  figure,  the  two  results  are  identical  because  the  three- 
dimensional  model  with  <p=  0  degenerates  to  the  two-dimensional  model.  The  dotted  line  in 

figure  3  shows  the  directivity  pattern  calculated  for  a  rigid  cylinder  with  an  equivalent  radius. 
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Figure  3.  Directivity  Patterns  at  f  =5000  Hz 
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To  compare  relative  scattered  pressure  fields  for  different  cases  (figure  3),  the  normalized 
peak  response  amplitude  scaling  factor  is  given  in  the  legend  for  each  case.  This  factor  was 
calculated  as  follows: 


Paf  =  201og10- 


<*>(  0,(p) 


9= 0 


fr  ) 

1/2 

j 

=  201og10<! 


xl/2 


nk0R^  cos  (p 


2X4B) 


n= 0 


dB. 


(101) 


Figures  4  through  12  present  directivity  patterns  calculated  for  the  coated  cylindrical  shell 
using  the  three-dimensional  model.  Figure  4  shows  the  directivity  patterns  calculated  at 
/=  5000  Hz  for  various  angles  of  incidence  <p.  In  figure  4,  the  solid,  dashed,  and  chain-dotted 
lines  denote  the  results  for  <p  =  0°,  30°,  and  60°,  respectively.  As  anticipated,  the  major  lobe 
width  of  the  directivity  pattern  becomes  broader  as  the  angle  of  incidence  increases.  In  the  limit 
the  directivity  pattern  becomes  omnidirectional  as  the  angle  of  incidence  approaches  90°. 

Figure  5  shows  the  directivity  patterns  calculated  for  the  angle  of  incidence  (p  =  45°  for 
various  frequencies.  The  solid,  dashed,  chain-dotted,  and  chain-dashed  lines  denote  the  results 
for/=  2500,  5000,  7500,  and  10,000  Hz,  respectively.  It  is  shown  in  this  figure  that  the  major 
lobe  width  of  the  directivity  pattern  becomes  broader  as  the  frequency  decreases.  In  the  limit 
the  directivity  pattern  becomes  omnidirectional  as  the  frequency  decreases. 

Figure  6  presents  the  effect  of  the  coating  thickness  h\  on  the  directivity  pattern.  The  solid, 
dashed,  chain-dotted,  and  chain-dashed  lines  denote  the  results  calculated  at /=  5000  Hz  and 
cp  =  45°  for  h\  -  2.54,  5.08,  7.62,  and  10.16  cm,  respectively.  It  is  shown  in  figure  6  that  the 
difference  between  the  results  is  not  significant,  although  the  results  show  some  different  values 
for  different  coating  thicknesses.  Figures  7  and  8  present  the  results  similar  to  those  shown  in 
figure  6.  Figures  7  and  8  show  the  directivity  patterns  for /  =  7500  Hz  and /  =  10,000  Hz, 
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respectively.  It  is  observed  in  both  figures  7  and  8  that  the  effects  of  the  coating  thicknesses  on 
the  directivity  patterns  are  not  significant. 

Figures  9  through  12  present  the  effects  of  the  coating  material  parameters  on  the 
directivity  patterns  calculated  for /=  5000  Hz  and  (p  =  45°.  Figure  9  shows  the  effect  of  the 
dilatational  wave  speed  (cdoi)  in  the  coating  on  the  directivity  pattern.  The  solid,  dashed,  and 
chain-dotted  lines  denote  the  calculated  results  for  the  dilatational  wave  speed  Cm  =  2000, 
20,000,  and  200,000  cm/s,  respectively.  As  shown  in  figure  9,  the  major  lobe  width  becomes 
narrower  as  the  dilatational  wave  speed  becomes  lower.  Figure  10  shows  the  effect  of  the  loss 
factor  (C,d\)  associated  with  the  dilatational  wave  speed  on  the  directivity  pattern.  The  solid, 
dashed,  and  chain-dotted  lines  denote  the  calculated  results  for  £/i  =  0.03,  0.3,  and  0.9, 
respectively.  As  shown  in  figure  10,  the  difference  between  the  results  obtained  for  different 
loss  factors  is  not  significant.  Figures  1 1  and  12  present  the  results  calculated  for  various  shear 
wave  speeds  (cj0i)  and  associated  loss  factors  (£h),  respectively.  In  figure  1 1,  the  solid,  dashed, 
and  chain-dotted  lines  denote  the  calculated  results  for  the  shear  wave  speed  cs0 1  =  2500,  5000, 
and  50,000  cm/s,  respectively.  In  figure  12,  the  solid,  dashed,  and  chain-dotted  lines  show  the 
calculated  results  £si  =  0.3,  0.9,  and  1.5,  respectively.  As  shown  in  both  figures  1 1  and  12,  no 
substantial  differences  are  noticed  for  various  shear  wave  speeds  and  the  associated  loss  factors, 
respectively. 
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Figure  4.  Effect  of  Incident  Angle  (p  on  the  Directivity  Pattern  atf=  5000  Hz 
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Figure  5.  Effect  of  Frequency  f  on  the  Directivity  Pattern  for  Incident  Angle  (p  =  45° 
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Figure  6 -  Effect  of  Coating  Thickness  hj  on  the  Directivity  Pattern  at  f-  5000  Hz  for  (p-45° 
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Figure  7.  Effect  of  Coating  Thickness  hi  on  the  Directivity  Pattern  at  f=  7500  Hz  for  cp- 45° 
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Figure  8.  Effect  of  Coating  Thickness  hi  on  the  Directivity  Pattern 
atf- 10,000  Hz  for  (p  =  45° 


Directivity  Pattern  atf  =  5000  Hz  for  (p  =  45° 
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Figure  10.  Effect  of  Loss  Factor  Qi  Associated  with  the  Dilatational  Wave  on  the 
Directivity  Pattern  at  f=  5000  Hz  for  cp  =45° 


Figure  11.  Effect  of  Shear  Wave  Speed  csoi  on  the  Directivity  Pattern 
atf-  5000  Hz  for  (p  =45  ° 
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Figure  12.  Effect  of  Loss  Factor  Qi  Associated  with  Shear  Wave  on  the 
Directivity  Pattern  atf=  5000  Hz  for  (p  =45  ° 
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CONCLUSIONS 


A  three-dimensional  analysis  has  been  made  of  the  scattering  from  a  coated  cylindrical 
shell  by  a  plane  acoustic  wave  making  an  incident  angle  with  the  normal  to  the  cylindrical  shell 
axis.  Based  on  the  limited  calculations  of  directivity  patterns,  the  following  conclusions  are 
drawn: 

1 .  The  major  lobe  width  of  the  directivity  pattern  becomes  broader  as  the  angle  of 
incidence  (with  the  normal  to  the  axial  direction)  increases  for  a  given  frequency. 

2.  The  major  lobe  width  of  the  directivity  pattern  becomes  broader  as  the  frequency 
decreases  for  a  given  angle  of  incidence. 

3.  The  major  lobe  width  of  the  directivity  pattern  becomes  narrower  as  the  thickness  of  the 
coating  increases;  however,  the  difference  between  various  thicknesses  of  the  coating  is 
insignificant. 

4.  The  major  lobe  width  becomes  narrower  as  the  dilatational  wave  speed  in  the  coating 
becomes  lower.  The  contribution  of  the  loss  factor  associated  with  the  dilatational  wave  to  the 
directivity  pattern  is  insignificant. 

Similar  conclusions  are  drawn  for  the  shear  waves  propagating  in  the  coating. 
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